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IMAGING SYSTEM 

FIELD OF THE INVENTION 

5 The present invention relates to an imaging system for body parts which utilises non- 
ionizing electromagnetic radiation, for example microwaves. In particular, although not 
exclusively, the imaging system is suitable for breast cancer screening. 

BACKGROUND TO THE INVENTION 

10 

Breast cancer is the most common cancer to affect women. The detection of 
malignancies at an early stage is deemed to offer the best prognosis for patients and this 
has lead to the establishment of screening programmes aimed at early detection. 

1 5 X-ray mammography is one commonly used breast cancer screening method due to its 
simplicity, high-resolution images and cost effective implementation. However, x-ray 
mammography has a number of associated limitations and drawbacks. X-rays are an 
example of ionizing electromagnetic radiation which can damage tissue and in some 
cases initiate malignant tumours. X-ray mammography requires the patient's breasts to 

20 be compressed between two plates which is uncomfortable for many women and makes 
it difficult to determine the true three-dimensional (3D) location of any suspicious 
features. Furthermore, women with silicone breast implants are also at risk from implant 
rupture due to the compression process. X-ray images are two-dimensional (2D) and a 
number of images from different views must typically be taken to provide some 

25 indication of the 3D location of suspicious features. X-ray detection of suspicious 
features relies on differences in density within the breast tissue under test and the 
density contrast between healthy and malignant breast tissue is small, typically only 
about 2%, which can make detection of tumours difficult. For post-menopausal women, 
x-ray mammography fails to detect up to 1 5% of cancers. For younger women, whose 

30 breast density is usually higher, up to 40% of cancers can be missed by x-ray 
mammography. Generally, the smallest tumour detectable with x-ray mammography is 
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about 4mm in diameter. A tumour this size is reckoned to have been in the body for 
about 6 years, that is, not particularly early in the tumour's development 

All of the above have provided significant incentive for researchers to develop 
5 alternative methods for breast cancer detection which alleviate some of the difficulties 
associated with x-ray mammography. Radar imaging, which utilises electromagnetic 
waves in the microwave region, has been identified as having potential for improved 
detection of breast cancer due to the large difference in complex permittivity between 
healthy and malignant breast tissue. US patent numbers 4,641,659, 5,807,257, 
10 5,829,437, 6,448,788, and 6,504,288 disclose various radar breast imaging systems. 

It is an object of the present invention to provide an improved imaging system for body 
parts, or to at least provide the public with a useful choice. 

1 5 SUMMARY OF THE INVENTION 

In one aspect, the present invention broadly consists in a method for generating a three- 
dimensional image of a body part including the steps of: 

scanning to obtain surface profile information for the body part; 
20 transmitting non-ionizing radiation through air toward the body part and then 

receiving non-ionizing radiation reflected back through air from the body part at a 
plurality of predetermined locations relative to the body part; 

obtaining radiation information at each of the predetermined locations from the 
reflected radiation received; and 
25 processing the surface profile information and radiation information to generate 

a three-dimensional image of the body part, the three-dimensional image having a 
plurality of image points. 

Preferably, the method further includes transmitting and receiving radiation at more 
30 than one frequency at each of the predetermined locations to obtain additional radiation 
information. More preferably, the method includes transmitting and receiving 
microwave radiation at a plurality of discrete frequencies between 10GHz and 18GHz. 
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Preferably, processing the surface profile information and radiation information 
includes synthetically focusing the radiation information to each image point of the 
three-dimensional image of the body part. More preferably, synthetically focusing the 
5 radiation information to each image point includes coherently adding radiation 
information from each of the predetermined locations. 

Preferably, the method includes moving an array of antenna elements relative to the 
body part and sequentially transmitting and receiving non-ionizing radiation from each 
10 of the antenna elements. Alternatively, the method includes sequentially transmitting 
and receiving non-ionizing radiation from antenna elements fixed at each of the 
predetermined locations. 

Preferably, obtaining radiation information at each of the predetermined locations 
1 5 includes measuring the reflection coefficient of the reflected radiation received. More 
preferably, this includes measuring the amplitude and phase of the reflection coefficient. 

Preferably, scanning to obtain surface profile information for the body part includes 
operating a three-dimensional laser profiler. More preferably, the method includes 
20 obtaining the surface profile information and radiation information simultaneously. 

Preferably, the method generates a three-dimensional radar image of the body part, 
which may be, for example, a human breast. 

25 In another aspect, the present invention broadly consists in an imaging system for 
generating a three-dimensional image of a body part including: 

a three-dimensional profiler arranged to scan the body part and obtain surface 
profile information; 

a radar device, displaced from the body part, arranged to transmit non-ionizing 
30 radiation through air toward the body part and then receive non-ionizing radiation 
reflected back through air from the body part at a plurality of predetermined locations 
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relative to the body part to thereby obtain radiation information at each of the 

* predetermined locations; and 

a control system arranged to operate the three-dimensional profiler and radar 
device, and which processes the surface profile information and radiation information to 
5 generate a three-dimensional image of the body part, the three-dimensional image 
having a plurality of image points. 

Preferably, the radar device includes a radiation source, radiation receiver, and one or 
more antenna elements which are operable to transmit radiation toward the body part 
1 0 and receive radiation reflected back from the body part. 

In the preferred form, the predetermined locations define a synthetic aperture over the 
body part within which the radar device obtains the radiation information. Preferably, 
the radar device includes an array of antenna elements moveable by an operable 

1 5 scanning mechanism, each antenna element of the array being selectively connectable to 
the radiation source and receiver via operation of a switching network. More preferably, 
the control system is arranged to operate the scanning mechanism and switching 
network to move the array within the synthetic aperture and sequentially operate each 
antenna element to obtain radiation information at each of the predetermined locations 

20 within the synthetic aperture. 

In an alternative form, the radar device includes a single antenna element moveable by 
an operable scanning mechanism, the antenna element being connected to the radiation 
source and receiver. Preferably, the control system is arranged to operate the scanning 
25 mechanism to move the antenna element to each predetermined location within the 
synthetic aperture for operation to obtain the radiation information. 

In another alternative form, the predetermined locations define a real aperture over the 
body part within which the radar device obtains the radiation information. Preferably, 
30 the radar device includes a plurality of antenna elements, one at each location within the 
fixed, real aperture. More preferably, the antenna elements are selectively connectable 
to the radiation source and receiver by a switching network and the control system 
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4 



operates the switching network to sequentially obtain radiation information at each 
predetermined location within the fixed, real aperture. 

Preferably, the radiation source is operable to generate non-ionizing radiation at more 
5 than one frequency. More preferably, the radar device is operable to obtain radiation 
information at more than one frequency at each predetermined location. In the preferred 
form system, the radiation source is a microwave source which generates microwaves at 
a plurality of discrete frequencies between 10GHz and 18GHz. 

10 Preferably, the radiation receiver is arranged to obtain radiation information by 
measuring the reflection coefficient of the reflected radiation received by the antenna 
element(s). More preferably, the radiation receiver is arranged to measure the amplitude 
and phase of the reflection coefficient. 

15 Preferably, the three-dimensional profiler is laser-based. Alternatively, the three- 
dimensional profiler is arranged to obtain surface profile information by utilising 
ultrasound or broadband microwave signals. 

Preferably, the control system is arranged to operate the three-dimensional profiler and 
20 radar device to obtain the surface profile information and radiation information 
simultaneously. 

Preferably, the control system includes an imaging system arranged to process the 
surface profile information and radiation information to generate the three-dimensional 
25 image of the body part. More preferably, the imaging system is arranged to synthetically 
focus the radiation information to each image point of the three-dimensional image of 
the body part. Even more preferably, the imaging system is arranged to synthetically 
focus the radiation information by coherent addition of the radiation information from 
each of the predetermined locations within the synthetic or real aperture. 

30 

Preferably, the control system is arranged to generate a three-dimensional radar image 
of the body part, which may be, for example, a human breast. 
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The invention may also be said broadly to consist in the parts, elements and features 
referred to or indicated in the specification of the application, individually or 
collectively, and any or all combinations of any two or more said parts, elements or 
5 features, and where specific integers are mentioned herein which have known 
equivalents in the art to which this invention relates, such known equivalents are 
deemed to be incorporated herein as if individually set forth. 

The invention consists in the foregoing and also envisages constructions of which the 
10 following gives examples only. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred forms of the invention will be described by way of example only and with 
1 5 reference to the drawings, in which: 

* 

Figure 1 is a rendered perspective view of a preferred form breast imaging system 
having a sensor head attached to a robot scanner; 

20 Figure 2 is a rendered perspective view of the sensor head of Figure 1 ; 

Figure 3 is a block diagram of the preferred form breast imaging system; 

Figure 4 is a block diagram of the radar device of the breast imaging system; and 

25 

Figure 5 is a schematic diagram showing the geometry relevant to the synthetic focusing 
algorithm. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

30 

The preferred form imaging system of the invention is a breast cancer screening tool 
and is arranged to scan a patient's breasts with microwave radiation in order to generate 
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3D radar images of each breast which can be examined for suspicious features such as 
malignant tumours. There is a large difference in complex permittivity between healthy 
and malignant breast tissue and this leads to greater scattered field amplitudes from 
malignant tumours embedded in healthy tissue which show up readily in a microwave 
image of scattered field intensity. For example, the real part of complex permittivity 
(the dielectric constant) for a malignant tumour is of the order of 50 at a frequency of 
10GHz whereas healthy tissue has a value of about 9. Hence, radar images are suited to 
breast tumour detection since the high permittivity contrast between malignant and 
healthy tissue translates to high-contrast images. 

The imaging system generates 3D radar images based on the intensity of the scattered 
field as a function of position from measurements of scattered fields external to the 
breast. The imaging system utilises a focusing algorithm to provide coherent addition of 
scattered fields at a given image point within the 3D radar image, thereby giving a 
measure of the scattered field intensity at a point in the breast being scanned. 

Referring to Figure 1, the preferred form imaging system 100 includes a sensor head 
101 which is translated relative to a patient 102 by a robot 103. The imaging system is 
arranged to scan each of the patient's exposed breasts individually and generate 
respective 3D radar images. In particular, the imaging system scans the patient's breasts 
to simultaneously obtain radiation information and surface profile information which 
are processed by an image generation algorithm to generate the 3D radar images. The 
sensor head 101 does not make contact with the patient 102 and there is no coupling 
medium, other than air, between the patient and sensor head during scanning. In an 
alternative form the sensor head 101 could be moved by means other than robot 103, or 
alternatively the patient could be moved relative to the sensor head. 

Referring to Figure 2, the sensor head is mounted to the robot scanning mechanism in 
the preferred form by a mounting flange 200. The sensor head includes a 3D profiler 
201 which is arranged to obtain geometric surface profile information of the breast 
during scanning. In the preferred form, the 3D profiler is a laser profiler device which 
uses a scanning laser stripe and charge-couple device (CCD) sensor to provide range 
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information by triangulation. The laser output power from the 3D profiler is deemed 
eye-safe. It will be appreciated that other types of 3D profiling devices could be utilised 
to obtain surface profile information of the breasts. For example, alternative forms of 
3D profilers may utilise ultrasound or broadband microwave signals to obtain the 
surface profile information. 

The sensor head also carries a radar device which is arranged to transmit non-ionizing 
radiation toward the breast and then receive radiation reflected back from the breast at a 
plurality of predetermined locations relative to the breast The radar device includes a 
radiation source 202 and receiver 203 which are connected to an array 204 of antenna 
elements or waveguides via a switching network 205. In the preferred form, the 
radiation source 202 is a Yttrium Iron Garnate (YIG) oscillator which generates 
microwaves over a broad range of frequencies and the radiation receiver 203 is a six- 
port reflectometer. The radar device is operated and controlled by an on-board computer 
system 206 and also has a calibration device 207 and an associated servo-motor 208. 

The preferred form radar device is arranged to obtain radiation information at a number 
of predetermined locations which define a synthetic aperture relative to the patient's 
breast. The radar device sweeps out the synthetic aperture by translating the antenna 
array 204 within the synthetic aperture and sequentially operating each of the individual 
antenna elements. For example, the preferred form radar device has a linear array of 
thirty two antenna elements arranged in two rows of sixteen antenna elements. During 
scanning, the antenna array is, for example, translated mechanically by the robot 
scanning mechanism to thirty two equally spaced locations in an orthogonal direction 
relative to the antenna array. At each of the thirty two locations, the thirty two 
individual antenna elements are sequentially connected to the radiation source and 
receiver by the switching network so that radiation information can be obtained at each 
of the 1024 locations of the synthetic aperture. In the preferred form, the antenna array 
has monostotic antenna elements, i.e. the antenna elements both transmit and receive 
microwave signals, but separate transmittal and receive elements could be used in an 
alternative arrangement. 
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The size of the synthetic aperture should preferably be no less than twice that of the 
body part to be imaged, so that the body part is illuminated sufficiently well by 
electromagnetic radiation from each antenna element. For imaging a breast, a value of 
15cm has been assumed as a typical linear dimension. Therefore, the minimum 
5 synthetic aperture size, D, is preferably twice this value, namely 30cm along each 
transverse axis. It will be appreciated that imaging system can alternatively operate with 
a smaller synthetic aperture to body part ratio depending on the system requirements. 

The required antenna element spacing in the antenna array is determined from the 
10 requirement to satisfy the Nyquist sampling criterion at the highest frequency of 
operation (shortest wavelength) so that grating lobes are avoided in the resulting image. 
This criterion requires that the element spacing be no greater than one half of a 
wavelength at the highest frequency of operation. For example, an upper frequency limit 
of 18GHz gives the largest allowed element spacing as 8.3mm. This element spacing in 
1 5 turn dictates the number of predetermined antenna locations in the synthetic aperture 
when combined with the minimum synthetic aperture size. 

The radiation information at each location within the synthetic aperture is obtained by 
illuminating the breast with microwave radiation from a transmitting antenna and then 

20 measuring the amplitude and phase of the reflected wave (scattered field) from the 
breast. The microwave radiation may be at a single frequency or small number of 
frequencies but in the preferred form radiation information is obtained at each location 
by repeating the measurement over a broad range of frequencies, one frequency at a 
time. In the preferred form radar device, a six-port reflectometer is incorporated into the 

25 microwave signal path. The six-port reflectometer is arranged to produce four voltages 
from diode detectors connected to its output ports from which it is possible to determine 
the amplitude and phase of the reflected signals relative to the incident (transmitted) 
signal. 

30 It will be appreciated that there are other alternative antenna arrangements which could 
be utilised to obtain the radiation information at each of the locations within the 
synthetic aperture. For example, the radar device may be equipped with only a single 
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antenna element which is translated mechanically to all locations with the synthetic 
aperture, although such an arrangement would be slow in terms of data acquisition 
speed. Alternatively, a real aperture could be provided in which there is an antenna 
element at each of the predetermined locations over the breast. With a fixed, real 
5 aperture the radiation information would be obtained by sequentially operating each 
antenna element one at a time. While a real aperture arrangement would be fast from a 
data acquisition viewpoint, it would also be more costly. The preferred form radar 
device is a combination of these two alternatives and is a compromise between data 
acquisition time and cost. 

10 

Referring to Figure 3, the sensor head 300 is mounted to a robot scanning mechanism 
301 and carries both the 3D profiler 302 and radar device 303. The robot scanning 
mechanism 301 is arranged to move the sensor head 300 relative to a patient's breast 
while the 3D profiler 302 and radar device 303 obtain surface profile information and 
15 radiation information respectfully as described above. A control system 304 is provided 
which controls the robot scanning mechanism 301, 3D profiler 302 and radar device 

303 during the breast scan. Further, the control system 304 is arranged to process the 
surface profile and radiation information with an imaging system to generate the 3D 
radar image of the breast. 

20 

Referring to Figure 4, the configuration and operation of the radar device 303 will be 
explained in more detail. The radar device 303 communicates with the control system 

304 via an on-board computer system 400. The radar device has a YIG oscillator 401 
which is operated in a swept frequency mode via its driver circuit 402 to generate 

25 microwave radiation at the desired discrete frequencies. The driver circuit 402 is in turn 
controlled by a sequence of binary signals from the on-board computer system 400. 

An important feature of the radar device is that the microwave power level emitted by 
each antenna element in the antenna array 403 is low and is of a non-ionizing nature. 
30 For example, the microwave power output from the YIG oscillator 401 may vary from 
30mW-50mW depending on the frequency. However, the power level made available to 
each radiating element in the antenna array 403 may be in the order of 0.1 mW due to 
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attenuation in the six-port reflectometer 404 and switching network 405. The sensor 
head 303 is also displaced, for example approximately 30cm, away from the patient's 
body which further reduces radiation exposure to the patient. Therefore, from a 
radiological stand point, the radar device is inherently safe. 

5 

Also, this stand-off distance is not critical but should preferably be greater than five 
wavelengths at the lowest frequency of operation in order that the illuminating 
wavefront from each antenna element has a spherical phase front with local plane-wave 
characteristics. That is, the breast is far removed from the reactive near-field region of 

1 0 the antenna and is illuminated by a wavefront having predictable phase and amplitude 
characteristics. A stand-off distance of ten wavelengths at the lowest frequency of 
operation is most preferable for reducing the effects of multiple reflections between 
breast and antenna which can contaminate the measured data and subsequent radar 
images. The stand-off distance is a compromise between being large enough to satisfy 

1 5 the above criteria and small enough that the transmitted and received signal levels are 
not too low due to the space-attenuation factor (that is the 1/R 4 dependence on the 
received power level, R being the object-antenna separation). This effect is 
compensated for in the preferred form by using a large number of elements in the 
synthetic aperture to enhance the received power levels when applying synthetic 

20 focussing. In addition, during the synthetic focussing process (which will be 
subsequently described), the size of the focal spot is also degraded (i.e. becomes larger) 
as the object-antenna separation is increased. To this end, it is desirable to maintain a 
focal ratio of the order of unity in determining the appropriate stand-off distance. 

25 A non-contact sensor head 300 enables the reflected signals from the breast to be 
accurately measured and allows calibration of the antenna system of the radar device in 
isolation. As mentioned, the stand-off distance between the breast and the plane of the 
synthetic aperture should preferably be at least 1 0 wavelengths at the lowest frequency 
of operation in order to reduce the effects of multiple reflections between antenna and 

30 breast to a negligible level. This allows the effects of the antenna system to be 
subtracted from the measured radiation information with the breast in place to give just 
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the reflectivity of the breast in isolation. A typical stand-off distance used for the breast 
imaging device is therefore 30cm at a minimum operating frequency of 10GHz. 

The radiation information to be measured by the radar device is the reflection 
5 coefficient of the reflected microwave signals at each location within the synthetic 
aperture and at each frequency of interest. In particular, the phase and amplitude of the 
reflection coefficient is measured. The six-port reflectometer 404 within the microwave 
signal path produces four voltages from diode detectors connected to its output ports 
from which the amplitude and phase of the reflected signals relative to the incident 
10 (transmitted) signal is determined. 

The six-port reflectometer 404 essentially combines the reflected microwave signal 
from the breast under test with a portion of the incident wave. This is done using four 
different relative phase differences introduced by the six-port reflectometer 404 between 

15 incident and reflected waves. The four combinations of microwave signals are then sent 
to four square-law detector diodes which generate four output voltages. One of the four 
output voltages is used as a reference such that three voltage ratios are derived for each 
measurement. These three ratios are converted into the real and imaginary parts of the 
reflection coefficient. The measured reflection coefficient information is then converted 

20 into digital data by an analogue-to-digital converter which in turn sends the digital data 
to the on-board computer system 400. 

The radar device employs a near-field imaging method in that the distance between the 
antenna elements and the patient's breast has a focal ratio typically in the order of unity. 
25 Therefore, the transmitted wavefronts illuminating the breast are highly curved. Further, 
the imaging system utilises an image generation algorithm which images objects 
embedded in the breast interior. In particular, the image generation algorithm takes into 
account the refraction at the various dielectric interfaces in order to focus effectively 
within the breast. 

30 

In the preferred form the radar device 303 includes a calibration device 407 and 
associated servo-motor 408 which are arranged to calibrate the six-port reflectometer 
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404 and antenna system. Calibration of the six-port reflectometer 404 will be described 
first. In order to accurately determine the complex reflection coefficient from the 
voltage outputs of the six-port reflectometer 404, it is necessary to calibrate the 
reflectometer to account for any imperfections and idiosyncrasies in the componentry. 
A number of 'calibration standards' are connected to the measurement port of the 
reflectometer and output voltages acquired as per a normal measurement. The 
calibration standards have known reflection coefficients for all frequencies of interest. 
For example, for the preferred form radar device, nine standards are used, all of them 
different lengths of short-circuited rectangular waveguide. 

It is possible to calibrate a six-port reflectometer using only five standards. However, a 
total of nine are made available in the preferred form imaging system due to the broad 
range of frequencies used. The key to an accurate calibration procedure for a six-port 
reflectometer is the selection of five standards with widely spaced reflection coefficient 
phase angles (the magnitude of the reflection coefficient is unity for all short-circuit 
standards). Having nine standards available allows one to select the five best phase 
angles for use at a given frequency thereby maintaining accurate calibration across the 
whole frequency band. 

The waveguide standards are built into the rotary calibration device 407 mounted on the 
sensor head which is able to connect each standard to the reflectometer measurement 
port one at a time by means of a servo-motor 408 to the appropriate position. 

For each of the nine calibration standards, the four six-port reflectometer output 
voltages are measured for each frequency in the band and stored. These are converted 
into real and imaginary parts of reflection coefficient and a set of calibration 
coefficients generated using a standard algorithm (not described here). The calibration 
coefficients characterise the six-port reflectometer 404 and enable the reflection 
coefficient of a breast under test to be accurately determined from the four diode 
detector output voltages taking into account the imperfections in the reflectometer 404 
itself. 
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The calibration of the antenna system will now be described. In order to extract the 
amplitude and phase of the reflection coefficients attributable purely to the patient's 
breast under test it is necessary to remove the contribution from the antenna system. 
This is done by performing a series of reflection coefficient measurements on the 
5 antenna system with no patient present. In particular, two measurements are carried out 
on the antenna system as outlined below. 

First measurement: With no patient present, the antenna system is positioned by the 
robot scanning mechanism so as to radiate into free-space with no reflective objects 

10 within close range. Each antenna element in the linear array is switched on in turn and 
the reflection coefficient determined for all frequencies via the output voltages from the 
six-port reflectometer. This represents the complex reflection coefficient of the antenna 
system and its associated switching network components and is referred to as the 
'empty room' case. The most significant contribution to the reflection coefficient in this 

15 case will be from the antenna apertures. 



Second measurement: The procedure outlined above is repeated with a metallic plate 
placed in close contact with the apertures of each antenna element in the linear array. 
This is referred to as the 'short circuit* case. The robot scanner moves the antenna array 
20 to a position where metal plate is automatically in close contact with the aperture 
plane. The most significant contribution to the reflection coefficient in this case will be 
from the short circuit plate. 

The short-circuit and empty-room data is used to extract the reflection coefficient of the 
25 breast alone from the overall measured reflection coefficient using the antenna array. 
This is an example of *de-embedding' applied to the measured reflection coefficient 
data to determine the reflection coefficient of the object in isolation. The de-embedding 
equation used is as follows: 



30 R bmul = R r* ...a) 



^empty Ix short 



where 
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Rbreast- Reflection coefficient of breast in isolation. 
Rmeas — Reflection coefficient of antenna and breast combined. 
Rempty - Reflection coefficient of isolated antenna radiating into free-space. 
5 Rshon " Reflection coefficient of antenna with short-circuit applied across apertures. 



The YIG oscillator 401 of the radar device generates continuous wave (CW) 
electromagnetic radiation covering a broad frequency bandwidth, i.e. the preferred form 
imaging system operates in broadband. In the preferred form, the operating frequency 

10 band is from 10GHz to 18GHz and radiation information is acquired at a number of 
frequencies throughout the band at each location within the synthetic aperture. The 
broadband frequency domain operation is utilised in order to provide a small focal spot 
size and hence good image resolution in the down-range direction. In the preferred form 
radar device, 161 discrete frequencies are used corresponding to a frequency interval of 

15 50MHz between 10GHz and 18GHz. The frequency interval is chosen to be small 
enough such that aliasing in the down-range direction is avoided in the final 3D radar 
images for the locations of interest in the image space. 

In order to obtain good focusing properties that approach the theoretical diffraction limit 
20 of half a wavelength for the size of the focal spot in the transverse plane, the synthetic 
aperture size needs to be large compared to the wavelength, X. Therefore, the 
requirement that D = 1 OA, at the lowest frequency (longest wavelength) follows. If D = 
30cm as mentioned previously, then X = 3cm. Therefore, the minimum frequency of 
operation for the imaging system is preferably 10GHz. 

25 

The broader the frequency bandwidth, the better the down-range resolution, so as broad 
a bandwidth as possible is desirable. Following on from above, 20GHz is preferably the 
upper frequency of operation. However, the vast majority of components will only work 
over a limited band, typically an octave at best. Therefore 18GHz is more practical from 
30 the performance of available components, giving a bandwidth of 8GHz. 
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The frequency interval between steps as the device is swept across the full frequency 
band is also determined by the need to satisfy the Nyquist sampling criterion. A small 
enough frequency interval needs to be used so as to avoid grating lobes in the time 
domain response resulting from an integration over the frequency domain data. This is 
5 in turn related to the round-trip time delay from source to receiver via the object under 
test. The frequency interval is chosen so that alias bands in the time domain response do 
not lie within the time interval for signals to make a round trip. This time delay can also 
be represented as an equivalent distance (there and back) in free-space referred to as the 
Alias-Free Range (AFR). A frequency interval of 50MHz is used in the preferred form 
10 breast imaging system giving 161 frequencies between 10GHz and 18GHz. 

Denoting the frequency interval by 5f, the corresponding separation of alias bands in the 
time domain, 8t, is given by the following equation: 



15 ^ = ~ .-(2) 



Equation (2) can be used to calculate an equivalent 'round-trip 9 distance in free-space, 
(AFR), by multiplying 5t by 2c where c is the speed of light in free-space to give 
equation: 



20 



AFR = 2cSt = — ...(3) 

df 



The microwave path length between source and image point and back should be less 
than the AFR in order to avoid contamination of the radar images from alias responses 
25 due to the sampling interval used in the frequency domain. Using 5f = 50MHz in 
equation (3) gives AFR = 1 1.99 meters in free space which is deemed to be sufficiently 
large for the proposed imaging system to avoid alias responses. A larger number of 
frequencies (and therefore a smaller frequency interval) could be used but this has to be 
offset against the total data acquisition time which must be kept small so as not to 
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inconvenience the patient. For example, the patient should ideally be able to hold there 
breath for the duration of the scan. 

The system described thus far gathers reflection coefficient data (radiation information) 
5 from the breast over a range of microwave frequencies. Images of scattered field 
intensity (3D radar images) are then generated by applying a synthetic focussing 
algorithm. 

Figure 5 shows the geometry of the antenna and breast configuration in a 3D Cartesian 
10 coordinate system. Referring to Figure 5 the vector Ri extends from the point in the 
antenna measurement plane to a point on the outer surface of the breast denoted by 
Ps(x s ,y s ,Zs). The vector R 2 extends from this point on the outer skin surface to a point on 
the interior skin surface. The vector R3 extends from this interior skin surface point to 
the image point P'(x\y',z'), the point at which microwave energy is to be focussed. 
15 This image point can be chosen arbitrarily. However, the path mapped out by the 
vectors Ri, R2 and R3 between antenna point and image point is not defined in an 
arbitrary fashion. Fermat's Principle is invoked so that the optical path is the minimum 
one possible. The minimum optical path, R m i n » is defined as follows for the geometry of 
Figure 5: 

20 

*nun = Minimum Value of {|R,| + V^T|R 2 | + V^«e |*M } »(4) 
where 

Sskin — Dielectric constant of skin. 
25 Stissue = Dielectric constant of breast tissue. 

There is one minimum path Rmi n for each image point and antenna point. So, for a given 
point in the image, there is a set of N Rmi n values where N is the number of antenna 
points used in the synthetic aperture. 

30 

The scattered electric field vector measured by the antenna at the point P(x,y,z) at a 
frequency denoted by the free-space propagation constant, k, is defined as i?scat(x,y,z,k). 
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The free-space propagation constant, k, is given by 2n/X where X is the free-space 
wavelength. A planar synthetic aperture is used here so that z = constant on the 
measurement plane. 

5 The radar image at a given point P* is now formed by applying a phase shift equal to 
2kRmin to the measured reflection coefficient data for each point in the synthetic 
aperture and then summing over all antenna locations. Summing over the frequency 
domain is also carried out. If the dielectric properties of the skin and breast tissue are 
assumed to vary negligibly with frequency (which is a good approximation), then the 
10 minimum paths between each image point and all antenna points will not depend on 
frequency. Therefore, once the minimum paths have been computed for a given 
combination of image point and antenna points, they can be used for all frequencies in 
the summation over the frequency domain. 

15 Mathematically, the above process can be represented by the following three-fold 
integral for generating the image, I, at P'(x\y\z'): 

/(*', /, z') = JJ *]E sca , (x, y, z, k)e 2JkR ^ dSdk ...(5) 

5 A, 

20 

where 

S = Synthetic aperture area. 

ki — Free-space propagation constant at lowest frequency. 
25 k2 — Free-space propagation constant at highest frequency. 

In equation (5) the factor of 2 in the phase shift term is present due to the need to 
account for the two-way path 'there and back' between antenna and image point. This 
phase shift term equalises the phase of the received signals from a given image point at 
30 all antenna locations so that when the summation over the synthetic aperture takes 
place, all quantities add up in phase to produce a much enhanced field at the image 
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point location. The measured fields are therefore focussed at the image point. This is an 
example of synthetic focussing applied to an antenna array. 

The use of the minimum optical path R min to calculate the appropriate phase shift is 
5 consistent with the Method of Stationary Phase often used to evaluate integrals of the 
type given in equation (5). This type of integral is characterised by a phase function in 
the integrand - often expressed as a complex exponential like that in (5) - which is a 
function of the integration variables. For values of the phase function which are varying 
rapidly with position, the oscillatory nature of the integrand in these regions results in a 

10 negligible contribution to the integral since positive and negative going portions of the 
oscillatory phase function tend to cancel each other out. The only significant 
contribution to the value of the integral comes from the region where the phase function 
is varying slowly such as in the vicinity of a stationary point in the phase function. This 
region corresponds to the minimum path Rmin and this is why it is used in the phase 

1 5 function exp(2jkR m i n ) of the integrand in (5). 

The vector nature of the electric field in (5) has been ignored since the dominant 
scattered field component will be co-polarised with the dominant polarisation present in 
the aperture of the antenna. That is, de-polarisation effects are ignored in the focussing 
20 algorithm - these will not be significant for a monostatic reflection coefficient 
measurement system. 

Equation (5) appears simple in form but complexity lies in the need to determine the 
values of Rmin for each combination of image point and antenna point. The 
25 determination of Rmin can be performed as a separate computational exercise and need 
only be computed once for a given antenna and breast geometry. In order to determine 
Rmin, it is necessary to have knowledge of the following: 

■ the geometric profile of the breast's outer surface relative to some known origin. 
30 ■ an estimate of the dielectric constants of the skin and interior breast tissue. 

■ an estimate of the skin thickness. 
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In the preferred form system, the geometric profile of the breast's outer surface is 
measured by the 3D laser profiler 201 co-mounted onto the radar sensor as described 
previously. Knowledge of skin thickness and dielectric constant of the skin and breast 
tissue to a high degree of accuracy is not necessary. An accepted value for the dielectric 
5 constant of skin at frequencies in the range 1 0GHz to 1 8GHz is 40 and that of the 
interior breast tissue is 9. The skin thickness may be nominally taken as 2mm. Values 
within 10% of the true values for dielectric constant will give rise to 5% errors in the 
optical path calculation due to the square-root dependence on the dielectric constants 
(see equation (4)). 

10 

For imaging purposes, the breast interior is assumed to be a homogeneous medium with 
a (mean) dielectric constant of s tisstie . While the breast interior will not be homogenous 
in practice, deviations from this mean dielectric constant will not be large for normal 
breast tissue. Large deviations from this 'background' dielectric constant - such as 

1 5 encountered with malignant tumours - will show up readily in the radar image whereas 
the smaller deviations in dielectric properties normally encountered with healthy breast 
tissue will scatter weakly and not show up as significant features in the radar image. 
Typically the imaging system of the invention will operate as a breast screening tool 
aimed at detecting the presence of suspicious objects within the breast rather than as a 

20 diagnostic tool. The above assumption of homogeneity for the breast interior is deemed 
sufficient for screening purposes. 

However, it is possible to glean additional information about the skin thickness and 
dielectric properties of the skin and breast interior from the focussed radar image data. 

25 The reflection coefficient due to the skin can be measured by focussing onto a region in 
the skin's vicinity and comparing the amplitude and phase of the measured result with 
that of a theoretical transmission line model. This is carried out at a single frequency. 
Data for multiple frequencies can then be generated in a similar manner and further used 
to determine the thickness of the skin, its dielectric constant and that of the breast 

30 interior (assumed homogeneous). Focussing at the skin surface helps to eliminate 
contributions from the chest wall giving values of reflection coefficient for the skin 
region only. 
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The minimum path Rmm is a function of the breast geometry as well as the antenna 
geometry and will therefore be unique to a particular patient. Values of R^in are 
calculated by fixing the antenna and image point locations and varying the position of 
5 the point P s on the skin's outer surface until the minimum value of the optical path is 
found. The two variables of interest here are x s and y s ,the x and y coordinates on the 
outer surface of the skin. The value of Zs is governed by the outer surface profile data 
(as measured by the laser system) and is a function of x s and y s . 

10 For a given point on the skin's outer surface, the point on the inner surface of the skin 
(where it meets the interior breast tissue) is automatically defined by Snell's Law of 
Refraction and so the vectors Ri,R2 and R3 are all fully defined for given values of 
antenna and image points along with values of x s and y s . Snell's Law of Refraction is 
wholly consistent with Fermat's Principle for a minimum optical path. Thus, the only 

1 5 variables in the search routine for the minimum path are x s and y s . 

Once found, the values of minimum path Rmm are stored in a five-dimensional array. 
Two indices are used to define the antenna location in the synthetic aperture and a 
further three to define the image point in 3D space. Image generation then proceeds by 
20 the numerical evaluation of the integral in equation (5). The image itself is usually 
displayed as the magnitude of the image function I(x\y%z'). 

Use of commercially available 3D visualisation software is the most effective means of 
displaying the 3D radar image data. Iso-surfaces and volume rendering visualisations 
25 are particularly appropriate for detecting suspicious features within the breast. 

The synthetic aperture method and apparatus described above consisting of an array of 
smaller antenna elements which behave collectively like an antenna of the same total 
physical size but whose characteristics can be reconfigured by manipulation of the 
30 relative phase and amplitude weighting applied to each element enables synthetic 
focusing to an arbitrary point in space via signal processing carried out after the data has 
been acquired in this piece-wise fashion. This provides a powerful microwave lens 
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which can be focused to an arbitrary location within the breast. This synthetic focusing 
ability provides the means of imaging small interior features such as malignant tumours. 
Also, due to the coherent addition of signals obtained from all elements in the synthetic 
array when focusing to a given point, the signal-to-noise ratio (SNR) of the 
measurement is improved by a factor N over a single measurement at a single frequency 
where N is the number of antenna elements in the synthetic array. Furthermore, by 
making measurements in the frequency domain, one frequency at a time, and then 
summing up the coherent signals from all antenna elements at all frequencies (to get a 
time domain response) the signal to noise ratio is further enhanced by a factor F where 
F is the number of discrete frequencies used. 

By coherent addition of signals at the designated synthetic focal point, the imaging 
device becomes very sensitive to scattered fields located at the focus. The coherent 
addition is carried out over all antenna locations and at all frequencies. A useful figure 
of merit is the increase in sensitivity of the imaging device as a result of focussing 
signals in this way and this is equal to the product of the number of antenna elements 
with the number of frequencies. This is also equal to the improvement in signal-to-noise 
ratio over and above a measurement of reflectivity carried out by a single antenna at a 
single frequency. For the breast imaging device this factor is 161 x 1024 = 164,864 
which is equivalent to an improvement of about +52dB. This is more than sufficient to 
overcome the two-way attenuation of signals in the breast tissue and skin which, at a 
depth of 5cm at a frequency of 18GHz, is about -40dB. To this end, higher frequencies 
than 18GHz could be contemplated with a subsequent improvement in resolution in 
transverse and down-range directions. 

In the preferred form frequencies in the range 10GHz to 18GHz are used. In general, 
attenuation in the breast tissue increases with increasing frequency. The benefit of using 
higher frequencies is the improved spatial resolution due to the reduced wavelength. 
The attenuation encountered does not pose difficulties for the preferred form method 
and apparatus of the invention due to the enhancement in sensitivity (e.g. +52dB) 
obtained as a result of coherent addition of received signals over a large number of 
antenna elements (e.g. 1024) along with integration over (e.g. 161) frequencies. Thus, it 
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is believed that the imaging system of the invention can accommodate higher 
microwave frequencies which enhances the resolution compared to lower-frequency 
systems. 

5 Also, the nature of electromagnetic scattering from objects small compared to the 
wavelength, such as the small malignant tumours of interest in breast cancer screening, 
needs to be considered. Such objects reflect incident energy back to the receiving 
antenna according to Rayleigh scattering theory. In Rayleigh scattering, the back- 
scattered power is proportional to the fourth power of the frequency. Therefore, the 
10 back-scattered signal from a small embedded object in the breast is 1.8 4 times larger at 
18 GHz than it is at 10GHz. This is a factor of approximately 10.5 or +10.2 dB. This 
enhanced scattering at the high-frequency end of the proposed frequency spectrum also 
helps to offset the increased attenuation in the breast tissue at the higher frequencies. 

1 5 The imaging system of the invention is also non-contact and does not require a liquid 
immersion medium surrounding the breast and antenna system. In addition, the 
separation between antennas and breast is typically of the order of ten wavelengths at 
the lowest frequency of operation (about 30cm at 10GHz). This is advantageous over 
some prior microwave systems which utilise both a liquid coupling medium and have 

20 antenna elements either in contact with the breast or in close proximity to it. The 
motivation for including a liquid medium around the breast is one of impedance 
matching with respect to the properties of the interior breast tissue. Reflections from the 
skin layer can be large thereby reducing the amount of energy entering the breast. If the 
dielectric constant of the liquid medium is similar to that of breast tissue then the 

25 amount of microwave energy penetrating the breast is maximised. The only residual 
effects that remain are reflections from the skin and attenuation in all media. 

Although the preferred form focussed synthetic aperture radar imaging system provides 
excellent sensitivity to the presence of internally embedded objects such as small 
30 tumours, further enhancement to the system's sensitivity is possible through the use of 
additional signal processing to reduce the contribution from skin reflections. The 
procedure used is to treat the breast interior as a homogeneous lossy medium for which 
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any reflections from the chest wall can be taken as much lower in magnitude than the 
reflections from the skin layer. This approximation is justified on the basis that the two- 
way attenuation of microwave energy in the breast tissue significantly reduces the 
magnitude of any chest-wall reflections as measured at the antenna. With this 
approximation in place, the breast interior can be regarded as infinitely deep in the z- 
direction and any reflections that take place are due entirely to the skin-tissue interface. 
The reflection coefficient from this idealised breast can then be estimated by applying 
the Physical Optics technique for electromagnetic scattering. The calculated scattered 
field is then subtracted from the measured data prior to applying the imaging algorithm. 
The effects of the skin are thereby significantly reduced by this procedure leaving only 
the contribution from the breast interior and any embedded features such as tumours. 

While the preferred form imaging system has been described as operating in the range 
of 10GHz-18GHz, the system could be arranged to operate within other higher or lower 
frequency ranges in the microwave band and the number of discrete frequencies utilised 
within the selected frequency range can be adjusted to suit design requirements. 
Furthermore, the imaging system may be arranged to operate at only a single microwave 
frequency, although this will degrade image quality. 

It will be appreciated that the aperture size within which radiation information is 
obtained can be altered as desired. Further, the number of predetermined measuring 
locations within the aperture and their respective spacings may be adjusted for specific 
requirements. For example, the number of predetermined measuring locations within 
the aperture may be increased to provide more radiation information in order to enhance 
the quality of the 3D radar image generated. 

While the preferred form imaging system has been described in the context of breast 
imaging, it will be appreciated that other body parts may also be imaged with the 
system, for example limbs to detect bone fractures or internal bleeding. Further, the 
imaging system of the invention could be adapted for brain scanning. 



210613-1 24 



Copy provided by USPTO from the IFW Image Database on 09/15/2005 



I 



It will be appreciated that the imaging system could be provided in the form of a hand- 
held portable scanning device which could be used in the field by ambulance drivers 
and the like. 

5 The foregoing description of the invention includes preferred forms thereof. 
Modifications may be made thereto without departing from the scope of the invention. 
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CLAIMS 

1 . A method for generating a three-dimensional image of a body part including the 
steps of: 

5 scanning to obtain surface profile information for the body part; 

transmitting non-ionizing radiation through air toward the body part and then 
receiving non-ionizing radiation reflected back through air from the body part at a 
plurality of predetermined locations relative to the body part; 

obtaining radiation information at each of the predetermined locations from the 
10 reflected radiation received; and 

processing the surface profile information and radiation information to generate 
a three-dimensional image of the body part, the three-dimensional image having a 
plurality of image points, 

15 2. An imaging system for generating a three-dimensional image of a body part 
including: 

a three-dimensional profiler arranged to scan the body part and obtain surface 
profile information; 

a radar device, displaced from the body part, arranged to transmit non-ionizing 
20 radiation through air toward the body part and then receive non-ionizing radiation 
reflected back through air from the body part at a plurality of predetermined locations 
relative to the body part to thereby obtain radiation information at each of the 
predetermined locations; and 

a control system arranged to operate the three-dimensional profiler and radar 
25 device, and which processes the surface profile information and radiation information to 
generate a three-dimensional image of the body part, the three-dimensional image 
having a plurality of image points. 
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